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ABSTRACT 


This  Investigation  studies  a  two-dimensional  model  of  a  layered 
medium  subjected  to  the  frictional  excitation  of  a  high-speed  asperity 
traversing  over  the  surface  of  the  medium.  The  general  analytical 
solutions  of  the  mechanical  stress  state,  the  temperature  field  and  the 
thermal  stress  state  are  obtained  and  expressed  in  the  Fourier 
transform  expressions.  Numerical  solutions  are  carried  out  for  cases 
of  uniform  and  parabolic  pressure  distributions.  The  resulting  stress 
state  yields  the  conditions  for  the  asperity  to  Initiate  cracks  In  the 
layered  medium.  This  paper  studies,  further,  the  thickness  effect  of 
the  surface  layer,  the  effect  of  relative  stiffness  of  the  surface 
layer  and  the  substrate,  and  the  effect  of  an  insulating  layer  versus  a 
conductive  layer. 
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CHAPTER  I 


INTRODUCTION 

1.1  Purpose  of  the  Investigation 

This  investigation  addresses  the  general  behavior  of  a  failure 
mechanism  that  is  caused  by  a  high  speed  asperity  traversing  over  a 
layered  surface  as  shown  in  Figure  1.  The  resulting  cracks  on  the 
surface  lead  eventually  to  total  failure  of  the  devices.  Such  a 
phenomenon  of  failure  has  been  observed  in  many  a  pair  of  mating 
surfaces  rubbing  against  each  other  such  as  brakes,  marine  seals,  and 
the  like. 

The  nominal  design  pressure  between  such  devices  is  based  upon 
the  total  mating  surface.  However,  at  the  operating  speed,  the  actual 
contact  area  can  be  smaller  by  a  factor  of  the  order  of  10-3,  or 
even  10~4.  Such  reduced  contact  area  could  be  the  result  of  an 
asperity  or  the  cause  thereof.  As  a  result,  a  low  design  pressure  may 
result  in  a  very  high  interfacial  pressure.  This  causes  a  large  fric¬ 
tional  force  in  the  contact  area.  The  high  friction  would  cause 
locally  an  extremely  high  temperature,  leading  to  cracking  of  the 
surface.  In  view  of  the  frequency  of  such  phenomenon,  many  attempts 
of  surface  modification  have  been  under  investigation.  Such  current 
significant  development  underlies  the  importance  of  the  present  inves¬ 
tigation  in  the  thermo-mechanical  behavior  of  layered  surface  under 
the  frictional  excitation  of  a  high  speed  traversing  asperity. 
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1.2  General  Background  and  Related  Investigation  in  Progress 


Two  bodies  in  sliding  contact  under  heavy  loads  will  experience 
high  local  temperatures  near  the  contacting  surfaces  because  of  exces¬ 
sive  frictional  heating.  Cracking  may  then  occur  in  the  neighborhood 
of  the  contact  zone  due  to  the  combination  of  thermal  heating  and 
mechanical  load.  This  phenomenon  is  known  as  heat  checking  [1], 

One  of  the  models  of  failure  is  due  to  high  pressure,  and  dry 
friction  exists  as  a  result  of  asperities  which  are  developed  on  the 
mating  surfaces.  The  asperity  may  be  either  an  external  or  internal 
material  inclusion  or  some  thermomechanical  deformation  of  the  mating 
surfaces.  The  asperity  represents  a  moving  concentration  of  pressure 
and  excess  Coulomb  friction  on  the  mating  surfaces.  Friction  influ¬ 
ences  both  the  mechanical  stress  state  and  the  thermal  stress  state 
because  frictional  heating  produces  a  severe  thermal  stress  state. 
These  high  stresses  can  then  initiate  a  fracture  of  the  material. 

A  general  survey  of  the  problem  of  cracking  through  the  develop¬ 
ment  of  a  frictional  hot  spot  was  discussed  by  Burton  [2],  The  exis¬ 
tence  of  high  hot  spots  was  demonstrated  experimentally  by  Sibley  and 
Allen  [3],  who  showed  systematically  moving  hot  patches  in  the  contact 
zone.  Two  and  three-dimensional  models  of  heat  checking  in  the 
contact  zone  of  a  bearing  seal  were  presented  by  Ju,  et  al .  [4], 

In  both  the  two-dimensional  and  three-dimensional  models  of  a 
single  moving  asperity,  the  high  temperature  field  is  concentrated  in 
a  thin  layer  In  the  neighborhood  of  the  surface.  The  depth  of  such  a 
surface  layer  is  of  the  order  of  twenty  percent  of  the  asperity  size. 
For  a  material  such  as  stellite  III,  the  high  temperature  field  will 
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not  alter  significantly  the  important  mechanial  properties.  Such  may 
not  be  the  case  for  other  materials.  When  property  changes  occur,  we 
will  have  a  problem  of  friction  cracking  involving  a  solid  with  a 
surface  layer  of  different  material  properties.  Problems  of  such 
phenomena  and  those  of  coated  surfaces  belong  to  the  general  class  of 
thermocmechanical  cracking  of  layered  media.  The  present  investiga¬ 
tion  deals  then  with  the  problem  of  heat  checking  in  media  with  a  hard 
surface  layer.  The  surface  layer  may  be  softer  or  harder  than  the 
substrate.  Their  thermal  properties  may  be  different.  These  factors 
and  the  thickness  of  the  surface  layer  constitute  the  parameters  of 
analysis  in  the  investigation. 

The  analytical  model  in  the  investigation  is  based  on  the  genei 
theory  of  a  continuum.  From  observation  of  failed  specimens  near  the 
surface,  it  is  reasonable  to  postulate  that  for  hard  surfaces  the 
plastic  deformation  and  rupture  at  the  surface  are  at  the  granular  or 
even  the  subgranular  level.  The  base  solid  material  subjected  to  the 
asperity  friction  is  essentially  elastic.  The  irrecoverable  work  in 
the  surface  deformation  manifests  as  heat  input. 

The  surface  layer,  be  it  thermally  induced  or  coated,  is  assumed 
to  be  of  uniform  thickness.  The  depth  of  the  substrate  is  infinite  in 
comparison  to  the  asperity  size. 

1.3  General  Theory 

1.3.1  The  uncoupled  theory  of  thermoelasticity  [5] 

The  basic  mathematical  formulation  of  thermoelasticity  describing 
the  behavior  of  continuous  media  are  the  following  equations: 
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(1) 


kT,ii  =  pcvT  +  (3x  +  2v)  aTo  ekk 
ai  j,  j  =  p‘U’i 

ei  j  “  7^ui ,j  +  uj,i ^ 

aij  =  5ijXekk  +  2ueij  '  5ij  (3x  +  2p)a(T  "  V  {4) 

The  coupled  dynamic  formulation  can  be  simplified  for  cases  without 
thermal  or  mechanical  shocks  to  the  uncoupled  quasi-static 
formulation. 

When  an  external  mechanical  agency  produces  variations  of  strain 
within  a  body,  the  time  rate  of  strain  variations  cause  a  change  in 
the  temperature  as  indicated  by  the  heat  conduction  equation.  Conse¬ 
quently,  by  a  flow  of  heat,  the  whole  process  increases  entropy  and 
therefore  increases  the  energy  stored  in  a  mechanically  irrecoverable 
manner.  This  phenomenon,  known  as  thermoelastic  dissipation,  is 
accounted  for  by  the  mechanical  term;  therefore  omitting  it  may  yield 
erroneous  results.  However,  the  dilatation  rate  in  equation  (1)  is  by 
order  of  magnitude  small  for  most  of  the  thermoelastic  problems.  The 
coupled  heat  equation,  [equation  {!)]  may  be  rewritten  as 


kT.ii  =  pcvT  1  + 


where  the  nondimensional  parameter  6  is  defined  by 
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(3X  +  2y)2a2To 


7T 


and  Ve  is  the  dilatational  wave  speed. 


(5) 


(6) 


5 


The  coupling  term  is  negligible  compared  to  unity  if 


For  illustration,  values  of  the  parameter  5  will  be  computed  for 
purposes  of  numerical  comparisons.  For  aluminum,  with 

X  =  4.13  x  1010  pa,  u  =  2.76  x  1010  pa, 

a  =  2.3  x  10~5/*C,  P  =  2.7  x  103  kgs/m3 
cy  =  835  N-m/(kg  #C), 

and,  taking  as  an  example  T0  =  90*C,  we  obtain  5  =  0.029. 

For  steel,  with 

X  =  12.4  x  1010  pa,  d  =  8.3  x  1010  pa,  a  58  1.2  x  10-5/°C, 

p  =  7.8  x  103  kgs/m3,  Cv  =  459  N-m/(kg  °C),  and  again  with 
T0  =  90°C,  the  corresponding  value  is  6  =  0.014. 

Thus  for  both  cases  the  coupling  will  be  small  [from  equation  (7)] 
if  approximately, 

£kk 

«  20. 

3aT 

For  temperature  distributions  with  no  sharp  variations  or  discon¬ 
tinuities  in  their  time  histories,  the  time  rate  change  of  the  dilata¬ 
tion  is  of  the  same  order  of  magnitude  as  that  of  the  temperature; 


thus  disregarding  the  coupling  term  as  described  previously  is 
reasonable. 

The  preceding  discussion  makes  it  clear  that  the  possibility  of 
omitting  the  coupling  terms  depends  not  only  on  the  fact  that  the 
inequality  6  «  1  must  hold  (as  it  does  for  most  metals),  but  also  on 
the  fact  that  strain  rates  must  be  at  most  of  the  same  order  of 
magnitude  as  temperature  rates.  The  latter  condition  implies  that  the 
time  history  of  the  displacements  closely  follows  that  of  the 
temperature;  in  other  words  no  pronounced  lag  or  vibrations  in  the 
motion  of  the  body  must  arise. 

1.3.2  Brittle  Fracture  [6] 

When  a  solid  is  subjected  to  increasing  loads,  the  resulting 
stresses  will,  at  a  certain  stage,  become  high  enough  to  cause  the 
solid  to  break  apart.  If  such  breakage  comes  about  before  the  piece 
has  thinned  down  to  zero  thickness,  it  is  called  fracture,  and  if  the 
amount  of  permanent  deformation  preceding  fracture  is  negligible,  it 
is  called  brittle  fracture.  In  this  investigation,  both  the  surface 
layer  and  the  substrate  materials  are  brittle  so  that  when  the  stress 
reaches  the  ultimate  stress,  brittle  fracture  will  occur.  Also,  if 
the  shear  stress  reaches  the  maximum  shear  stress  at  the  interface, 
shear  delamination  will  happen. 
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CHAPTER  II 


MATHEMATICAL  MODEL 

The  problem  under  consideration  is  one  of  thermomechanical  crack¬ 
ing  in  a  layered  medium  by  a  fast  moving  asperity,  as  shown  in  Figure 
1.  Because  the  actual  contact  area  of  the  moving  asperity  is  much 
smaller  than  the  layered  medium,  the  mathematical  model  is  represented 
by  a  half-space  subjected  to  a  fast  moving  asperity  whose  effect  is 
delineated  into  a  moving  heat  source  and  a  moving  mechanical  load  of 
combined  pressure  and  tangential  friction.  Two  sets  of  coordinates 

I  I 

are  employed:  x^  -  x2  are  fixed  to  the  medium,  xY  -  x2  are  fixed  to 
the  moving  load.  The  relative  speed  of  the  contact  surface  is  assumed 
to  be  large  enough  to  result  in  a  high  Pdclet  number  (vt/ie)  but 
smaller  than  the  Rayleigh  wave  speed.  The  quasi-static  theory  thus 
holds.  The  general  solutions  for  the  mechanical  stress  state,  the 
temperature  field  and  the  thermal  stress  state  are  approached  by  the 
use  of  the  Fourier  transform  method.  Because  of  their  complexity,  the 
general  solutions  are  left  in  the  transformed  space.  Solutions  by 
numerical  integrations  of  the  inverse  transforms  are  carried  out  for 
special  problems. 

2.1  Basic  equations 

The  governing  equations  come  from  Cauchy's  law  and  the  uncoupled 
theory  of  thermoelasticity,  in  terms  of  the  moving  convective 
coordinates  { } .  The  acceleration  in  Cauchy's  law  will  have  only 
the  convective  terms.  Hence,  Navier's  equation  is  expressed  as 
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(x  +  u)  uk  kl-  +  *  (3X  +  2y)  a  T,.  =  pV 


2  3  ui 


3X 


1 


where  summation  convention  is  used  for  repeated  indices  of  roman 
minuscules. 

For  i  =  1 

2  2  2 
«  3  U.  3  1*2  3  U. 

(X  +  2u  "  pv  )  — y  +  (X  +  u)  ■  ■  +  u  — 3 (3X  +  2y)  a  -re—  *  0  . 

3xJ  3X13X2  3X|  3X1 

For  i  =  2 

2  2  2 

p  3  U?  3  U.  3  Up  y 

(u  “  pv  )  — y  +  (X  +  p)  gx  3x—  +  (X  +  2y)  — ^ - +  2y)  a  jy  -  0 

3x^  1  2  3Xp  2 

The  Hooke’s  law  equations  are 


3u.  3Up 

an  =  (x  +  2u)  +  x  *rz "  {3x  +  2v)  a(T  "  To)  • 


3Uj  3Up 

a12  "  v  \ 3X^  +  Tx^  /  * 


3u,  3Up 

°22  =  X  3x7  +  (X  +  2m)  3^  ‘  {3X  +  2w)  °{T  "  V  * 

The  above  equations  apply  to  both  the  surface  layer  and  the  substrate. 
The  surface  boundary  is  traction  prescribed  so  that 


'12^  3  wf  p(xl^  (  xi  i  *  »  x2  =  °)  » 
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u> 

J22 


It 

o 

( 

xll 

(-P(Xj) 

( 

xl| 

>  *  »  *2 
<  l  ,  x, 


0) 

0) 


Regularity  conditions  hold  at  infinity. 


°11*®12,022*U1*U2  +  0 

(°  £  x2  —  H 

® \2. 9  ^22  *  ^  1  *  ^2  **  ^ 

0  p 

(x^  +  x2  + 

) 


«) 


Continuity  conditions  apply  at  the  surface  layer/ substrate  interface 


(1L 
12  “  °12  * 


H) 


where  the  superscript  (1)  designates  the  surface  layer,  (2)  the 
substrate. 

The  heat  conduction  equations  are 


S2T(1)  ,2tU) 

1  ~~i 


ax. 


ax. 


aT 


(l) 


3xi 


32t(2)  32t(2) 
T”  ~TT~ 


ax. 


ax; 


_  v  aT 


(2) 


k2  3xl 


The  boundary  conditions,  the  continuity  conditions,  and  the  regularity 
conditions  for  the  heat  equations  are 

(  <.  *  ,  x2  *  0) 

<|xl|  >  ‘  •  *2  '  0) 


-k 


aT 


(l) 


l  ax. 


aT(1)  =  .  3T(2) 

ax2  ~  k2  ax2 


H) 
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T(l)  =  t(2) 


(x2  =  H) 


T<»  .  I(2) 


(x2  +  x|  ♦  «)  . 


The  following  dimensionless  quantities  will  be  used  for  general 
computation. 

5  ~  Xj/ 1  ,  n  _  x^/ 1  i  D  *  H/t  i  -  oj^/Pq  »  0^  ~  a2Z^o  * 


*(i)  =  (T(i)  -  T  )  k1/qQt  »  Q  *  g~  »  Rx  =  vt/^  ,  and  R2  =  vt/^  . 

The  linear  theory  of  elasticity  allows  the  equation  to  be  delineated 
into  the  mechanical  and  the  thermal  parts.  They  are  thus  grouped  into 
(1)  equations  for  the  mechanical  stress  field,  (2)  equations  for  the 
temperature  field,  and  (3)  equations  for  the  thermal  stress  field. 
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2.1.1  Mechanical  stress  field 
In  the  surface  layer 


(2)  _  w2 


nn 


r* 


or  -  21 


3U 


(2) 


35 


+  v 


2  3V 


(2) 


3n 


(17) 


The  boundary  conditions,  the  regularity  conditions,  and  the  continuity 
conditions  are  correspondingly 


=  wf  P(^)  (| 5  |  1  1  »  n  =  0)  , 


(1). 

j 

nn 


(5> 


(|e|  >  1  .  n  -  o) 
(|e|  i  1  .  n  -  0) 


(18) 


(19) 


55 


,o„  ,a  ,u,V  +  0  (0<n<D>  H  +  »), 

5n  nn  -  —  I 


(20) 


arF»0F««° nn*U*V  +  0  ^  +  *  *  *)  » 

55  5n  nn 


(21) 


(1)  .  o‘2),  o(1)  -  „(2),  a<u  •  U(2),  V(1)  .  V(2>  {„  -  D)  .  (22) 


'En  En  nn  nn 


2.1.2  Temperature  field 

(.V11  ,  a24(1)  .  p  an111 

IT  ^  ’ 


(23) 


££L  +  jV!i  .  „  a±!!i 

as2  17“  2  35 


(24) 
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The  boundary  conditions,  the  continuity  conditions,  and  the  regularity 
conditions  are 


36U)  *  ' 

|Q(€)  (|c  l  <  1  .  n  =  0) 

(o  (|c  |  >  1  ,  n  =  0) 

(25) 

n!fi 

3n  w  3n 

(n  -  D)  , 

(26) 

♦“>  *  *(2) 

(n  -  D)  , 

(27) 

♦u).*<2)  *0 

U2  +  n2  -  -)  . 

(28) 

2.1.3.  Thermal  stress  field 
In  the  surface  layer 


*2 


Y1  3± 


(1) 


,.2  .2.  32u(1)  / .2  ,2,  32v(1)  ,2  32u(1)  b 

“IT”  '  T  +  J  —2 - 

n  c2  (29) 


«  0  , 


*2 


9  9  *2„(1>  9  9  *2(1)  7  _ 2V ( 1 )  b 

(J2  -  M2)  ?  %■  -  +  (I2  -  J2)  3-“  •  +  I2  3  V 

3€2  H3n  3n'  c 


(1) 


Y1  34  -  n 

r~  15  3n  ~ 0  • 


(30) 


5po 


,2  3U 


(1) 


H 


+  (I‘ 


2J2)  3V 


(1) 


an 


*2 

Y1  JD 

~r  * 


(31) 


,(1) 

Cn 


2  u2 

«02  y- 

o 


(1)  ,u(l) 


3U  ,  3V 


3n  3C 


(32) 


14 


nn 


0 


(I2  -  2J2)  3u 


(1) 


n 


+  i 


2  3V 


(1) 


h*2 

5  Y1  .(1) 


3n 


(33) 


In  the  substrate  region 


2  2  a2u(2)  2  a2V(2)  32U(2)  b  y2  3*( 2 } 

(N^  -  M2)  3  UA  +  (N2  -  1)  — +  -2JL- - -i  14 = 

1?^  353*  ^7^  If  35 


0  ,  (34) 


2  a2V  2  2  32u(2)  2  32V(2)  b  y2  3d,(2) 

(1  -  M2)  ■  %—  +  Ur  -  1)  +  N2  - 14 -  =  0  ,  (35) 

n  35311  3n2  C2  3n 


(2) 


U2 

P. 


2  an ^ 2 ^  2  3V^2^  b  y2  (2) 

N2  +  (N2  -  2)  3V  -  2 


3C 


3n 


(36) 


(2)  _  v2  f3u(2)  3V(2)~| 

°£n  "  P0  |_  3n  3C  J  ’ 


(37) 


(2)  _  w2 


nn 


(N2  -  2)  9U 


(2) 


*  N 

35  3n 


2  3V(2)  b  y2  (2) 


(38) 


The  boundary  conditions,  the  continuity  conditions,  and  the  regularity 
conditions  are 


«  o 

5n 

(n  =  0)  , 

(39) 

»  0 

(n  *  0)  , 

(40) 

nn 

5C,aCn,ann’ 

u,V  0  (0  <  n  <  D  ,  5  -*•  •)  , 

(41) 
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,  r2  .  2 
U  +  n 


■) 


(42) 


<W‘VU>V  *  0 


(1)  _  (2)  (1)  .  (2)  u(l)  _  (2)  (1)  .  (2)  ,  _  ni 

aZn  ‘  %  *  %  '  ann  *  u  "  u  *  v  "  v  (n  -  D)  . 


(43) 


2.2  General  solution 


The  Fourier  transform  [7] 


F{  }  =  (2*)_1/2  /  {  }  eisc  d5 

”*oo 

is  used,  denoting  the  transformed  quantities  by  a  superposed  tilda  -  , 

then  letting  u  =  iu,  v  =  v,  ;  =  Sn,  »  a.J s,  or  =  1op  /s,  a  * 

CE  CC  En  Sn  nn 

a  /s,  P  =  P/s,  i  =  i/s,  Q*  =  Q*/s,  and  1  =  4-  =  -  4—.  We  can  obtain 
nn  a?;  s  an 

the  analytical  solutions  in  the  transformed  space. 


2.2.1  The  analytical  solutions  of  mechanical  stress  field 
Equations  (8)  through  (17)  become 


02  J<‘>"  -  (I2  -  M2)  S'11  * 

(I2  -  J2)  V(1)'  =  0  , 

(44) 

I2  V*11  -  (02  -  M2)  v(1)  - 

(I2  -  J2)  u(1)’  -  0  . 

(45) 

;‘J>  -  5^[-i2  :UI  ♦  (I2 

-  2J2 )  V{1)]  , 

(46) 

;'l»  ■  $j2  %  (utu'  +  v(u) 
5n  po 

> 

(47) 

16 


(48) 


;(1)  =  «£  n2  v(1)‘ -  (i2  -  2J2) -u(1)]  , 
nn  P0 

(48) 

u(2)“  .  (n2  _  m2j  -(2)  +  (fJ2  _  y( 2) '  =  0  j 

(49) 

N2  V{2)"  -  (1  -  M2)  V{2)  -  (N2  -  1)  u(2)'  =  0  , 

(50) 

Op2)=  %  [-N2  u(2)  +  (N2  -  2)  V(2)‘]  , 

(51) 

(u(2)‘  +  v(2))  , 

(52) 

-(2)  =^2[n2  y(2) 1  _  (fJ2  .  2)  -(2)]  # 

nn  P. 

(53) 

The  boundary  conditions,  the  regularity  conditions  and  the  continu¬ 
ity  conditions  are  also  in  transformed  expressions 

i  r 

=  iufP  -  f  P*(C)  e1s?  d?  at  c  =  0  (54) 

T  /2ns  J 


o(1)  =  -  P  = - L_ 

nn  m  s 

where  P*(?)  =  | P  ^ 

I  0 


P*(5)  e1SC  d£  at  5  =  0 


(55) 
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°SE’  °W  U>  V  *  °  (‘ 


(56) 


-  (1)  -  (2)  -  (1)  -  (2)  -(1)  =  -(2)  w( 1)  =  y(2) 

°5n  °5n  ’  ann  %  ’  U  U  ’  V  V  (57) 

(c  =  SD) 

From  equations  (49)  and  (50)  we  can  get 
ul2)  =  f‘2  (f2  -  g2)""1  [ (1  -  f2)  V(2),,‘  +  (f4  -  g2)V{2)']  (58) 


5(2) 


(4) 


(f2  +  g2)  V{2)"  ♦  fV  V(2}  »  0 


V'-'  -  If-  +  g“j  v  +  rg" 

where 

f  =  (1  -  M2/N2)1/2  ,  g  =  (1  -  M2)172  . 


(59) 


For  the  subsonic  case 


M  =  v/c‘  <  1,  N  =  c1/c2  >  1 

,  9  1/2  _  „  1/2 

f  =  (1  -  tr/ir)  =  (1  -  vVCj  )  >  0 

g  -  (1  -  mm  -  (1  -  v2/c22)1/2  >  0 


and 


Cj  >  c2  so  that  f  >  g. 


Equations  (44)  and  (45)  yield 


;(l>  -  u2  -  m2)'1  (i2  -  o2)"1  i2j2v(1)"' 


+  [(I2  -  J2)  -  J2(J2  -  M2)]  V(1)l 


and 


v(1)(  -  C(i  -  m2/j2)  +  (i  -  m2/i2)]v(1)" 

+  (1  -  M2/!2) (1  -  M2/J2)Vll)  =  0 


(60) 


(61) 
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Equation  (61)  has  a  characteristic  equation 


X4  -  [(1  -  M2/J2)  +  (1  -  M2/ 12 ) 3 X2  +(1  -  M2/I2)(l  -  M2 / J 2 )  =  0 


in  which  x  is  the  characteristic  value  of  =  ex?. 
For  the  subsonic  case 


M2/J2  <  1,  M2/I2  <  1,  and  I  >  J. 

with  j  =  (1  -  M 2/I2)1^2,  k  =  (1  -  M2/J2)1{2  then  applying  the  boundary 

conditions,  we  can  obtain  the  solutions  for  s  >  0. 

-(1)  *  *  *  * 

V  -  ax  cosh  jns  +  bj^  sinh  jns  +  c1  cosh  kns  +  d  sinh  kns  ,  (63) 


,-i  *  *  *  ★ 

j  (ax  sinh  Jns  +  bx  cosh  jns)  +  k( Cj  sinh  kns  +  dx  cosh  kns). 


y  [~6(M2j  1  +  2J2 j)  (a1  sinh  jns  +  bj^  cosh  j^s) 


-  25J2k(c1  sinh  k^s  +  d*  cosh  kps)]  ,  (65) 


.(in  Up  2  *  *  9  9 

o'  =  y  [25J  Uj  cosh  jns  +  b1  sinh  jns)  +  5J  (1  +  k  ) 


(Cj^  cosh  kns  +  dj^  sinh  kgs)]  ,  (66) 


.[  1  \  Vo  O  .1  0  it  ★ 

o'  a  3=-  [SJ  j  (1  +  k*1)  (a,  sinh  jns  +  b.  cosh  j^s) 
nn  l  l 


C  I 


(68) 


V(2)  =  A*e“fn S  ♦  B*  e‘^S 


(2)  =  -  gB;  e-^s  . 


(69) 


Z{2)  =  y-  Cf_1(M2  +  2f2)  A*e"fnS  +  2gB*  e"9nS]  , 


(70) 


;<«-  £  C2A;  e~fns  ♦  (1  *  g2)  sj  e'”"5]  , 


(71) 


;(2)  -  %  C-f'1!!  *  g2)  A?  e'f”s  -  2g8;  . 

nn  rn  i  l 


(72) 


★  ★  ★  ★ 

For  s  <  0  we  can  get  a  similar  set  of  solutions.  Where  alt  blt  c1#  dx 

★  ★  .... 

Aj,  and  B2  are  functions  of  s  and  depend  on  the  pressure  distribution 

profile  P,  that 


P  = 


2  sin  s 


✓2*  s 


T 


for,  uniform  pressure  and 


P  =  (S2LJL  -  cos  s) 

/2w  sJ  s 


for  parabolic  pressure. 
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Now  enforcing  the  boundary  conditions  (54),  (55),  and  (57),  the 
★  ★  ★  ★  ★  ★ 

six  constants  A^  alf  blt  c^  and  dt  are  interrelated  by  six 
algebraic  equations 


*  ,  2  *  - 
2h  *  «  +  k  >  ci =  ^7  2 


-1  ,  9  *  *  P 

j  1  (1  +  k*)  bx  +  2kd:  =  -  -I-? 

5  J 


-1  *  -1  ★  *  *  -I  -fsD  * 

J  Sj  ix  +  j  c.  b,  +  kstc1  +  kc^d,  +  f  e  A,  +  ge 


'j  U1  k  1 


*  *  *  *  -fsD  *  -asD  * 

cjal  +  sjbl  +  ckcl  +  skdl  ‘  e  A1  '  e  B!  =  0 


*  *  7  *  7  *  7b. 

2cja1  +  2Sjb1  +  (1  +  k^)ckc1  +  (1  +  r)s0d 


■fsD 


k  1 


5J 


(1  +  q2)e"9S°  * 

-  - 2-IH -  B.  =  0 

SJ*  1 


-1,  2.  *  -1-  2  *  *  * 

j  (1  +  k  )Sja1  +  j  (1  +  k  JCjbj^  +  2ks|cc1  +  2kck  dj. 


♦  U  *  92)  e-fS0  .*  *  i2i-’SD  -* 


fSJ 


50 


Bi =  0 


(73) 

(74) 


"9SDB*=  0  (75) 
(76) 


* 


(77) 


(78) 


In  the  above 
s.  -  sinh  (jsD) 

J 

c.  =  cosh  (jsD) 

J 

s^  =  sinh  (ksO) 
c^  =  cosh  (ksD) 


represent  the  six  algebraic  equations  in  matrix  form,  then  solve  it  by 
Cramer's  rule. 
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The  expressions  of  a,  A}»  a2,  A3,  A4 ,  A5,  and  a6  are  given  in 
Appendix  I. 


22 


2.2.2  The  analytical  solutions  of  temperature  field 


Equations  (23)  and  (24)  in  the  transformed  space  became 
;(1)"  -  (1  -  TRj/s)  *(1)  =  0  ,  (79) 

*(2)”  -  (1  -  iR2/s)  *(2)  =  0  .  (80) 

Applying  boundary  conditions  we  can  get  the  solutions  for  s  >  0. 

m  i  -/I  -  iR,/s  ns  «  / 1  -  ift./s  ns 

■  Ax  e  L  +  Bx  e  1  ,  (81) 


I 


-/I  -  iR2/s  ns 


(82) 


I  I  I 

where  Alf  Blf  and  A2  are  functions  of  s  and  depend  on  the  heat  input 
profile  Q*,  that 


_  2  sin  s 

y - 5— 

m  S2 


for  the  uniform  pressure  and 

Q*  3  — - — T  (il£-i  -  cos  s) 
m  s3  3 

for  parabolic  pressure. 

For  s  <  0,  Equations  (81)  and  (82)  may  be  used,  provided  s  is 

-12) 

replaced  by  n  =  -s  and  a  negative  sign  is  added  to  y  . 

Apply  the  boundary  conditions  and  the  continuity  condition  at  the 

I  I  I 

surface  layer/ substrate  interface,  At,  Bj,  and  A2  can  be  readily 
solved. 
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*:  -£ 


F  +  bG 


1  S  F2(l  -  e-2FsD)  ♦  BFGd  ♦  e‘2FsD) 

’  u  £  (F  -  bG)  e~2FsD 

1  s  F2(l  -  e'^FsD)  +  sFG ( 1  +  e“2FsD) 

■  _  Q*  2e(-F+G)sD 

2  s  F(1  -  e~^s0)  +  bG(1  +  e"2FsD) 


where  F  =  (1  -  1R1/s)1/2  ,  G  =  (1  -  iR2/s)1/2 


2.2.3  The  analytical  solutions  of  thermal  stress  field 
Equations  (29)  through  (38)  become 


J2u(l)"-  d2-M2)  ;(1)  +  (I2-  J2)  v(1)'  . 


I2V(1)"  -  (J2  -  M2)  V(1)  -  (I2  -  J2)  u(1)'  =  — i  *(1)'  , 


a(1)  =  5 


u2  2  -(1) 
y  -r  UU'  + 

0 


(I2  -  2J2)  V(1)'  -  *(1)  , 

C2 


:»>.«* 5 (!»»'♦  ¥»>) . 


;(1)  .  5“2  ,25(2)'  .  (I2.  2J2,  ;<1)  .fin  ;<»>  , 

no  r  c 


(88) 


u(2)'  -  (N2  -  M2)  U(2)  +  (N2  -  1)  V(2)'  =  —2-  $(2)  , 


iZ  V(2)"  -  (1  -  M2)  V(2)  -  (N2  -  1)  u-{2)'  i(2)‘  , 


Op2)  =  -N2  u(2)  +  (N2  -  2)  V(2) 


5  y2  -(2) 


;(2L^(-(2)'  +  v-(2)) 
*\> 


-(2) _  ^2 
*nn  P„ 


-  (N2  -  2)  ul2)  +  N2  V{2)'  - $(2) 

C2 


The  boundary  conditions,  the  regularity  conditions  and  the 
continuity  conditions  are 


o  «t  =  o 


o(1)=  0  at  c  *  0 
no 


aP,i  o_  ,  o  ,  u,  v+0  when  (95) 

5n  nn 

-(1)  _  -(2)  -(1)  _  -(2)  -(1)  _  -(2)  .  s(l)  _  j(2)  _  _  - 

0„  =  0,.  ,  o  =0  ,u  =u  ,  and  V  =  V  at  c  =  sD 

Sn  5n  nn  nn 

-(1)  -(1)'  —(2)  -(2)' 

where  ♦  ,  <)>  ,  4>  and  $  come  from  temperature  field. 


;<!>•  ,  F(.A;  e-Ft  t  B;  eFC  , 
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;l2,‘  -  -  G2a;  e”Gc 


A^,  Bif  and  A£  are  the  same  as  l'n  the  temperature  field.  Equations 
(88)  and  (89)  yield 

u(2)  =  f“2  (f2  -  gV1  [(1  -  f2)  V(2),,‘  +  (f4  -  g2)  V<2)‘]  (9 

-  E^(2)  -  e2*(2)" 


f(2)  .  (f2  +  g2,  y ( 2 )  +  f2  2  y(2)  ..  -(2)* 


v-  +  r<r  v-  -  e3»-  -  e4  (98) 


1  u  Tp  7  D  Y; 

where  E.  =  — *-= — » - *=■  ,  E,  =  — * - if — * - , 

1  (N2  -  M2)  c2  (N2  -  M2) (N2  -  1)  c2* 

2  h2  K2 

E  =  T2  ,  £  -  1  Y2 


4  7  c. 


;l2)"  -  g2a!  e’G< 


and  i(2)"'  =  -G3A2  e'Gc 


Equation  (98)  is  a  fourth  order  nonhomogeneous  ordinary  differential 
equation,  where  the  nonhomogeneous  part  comes  from  temperature  field. 
The  complementary  solution  of  equation  (98)  is 


V<2>  =  C  +  g*gf?  +  C2  +  D^e^1* 


Assume  the  particular  solution  of  equation  (98)  is 


V(2)  =  A^ 


(100) 
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Substituting  (100)  Into  (98),  we  can  obtain  A|  . 


a!  -I- 


(-E3G  +  E4Gj) 


2e 


(-F+6) sD 


3  S  [G4  -  (f2  +  g2)G2  +  f2g2]  [F(l  -  e~2FsD)  +  eG(l  +  e"2FsD)] 


The  solution  of  equation  (98)  is  the  combination  of  the  complementary 
solution  and  the  particular  solution 

V(2)  =  A*e"fc  +  B*e‘9?  +  c2ef5  +  D^e95  +  A3e"G?  (101) 

Equation  (101)  together  with  the  regularity  condition  (95)  yields  a 
solution  for  s  >  0. 


3«)  „ 

2 


V'“'  =  A„e'f?  +  B*e‘9c  +  A„e 


*  -Gc 
3" 


(102) 


Substituting  (102)  into  (97)  u^  can  be  obtained. 


u{2)  *  -f”1  A*e“f?  -  gB*e“95  +  Ece"Gc 


(103) 


where  Es  =  -  \  f ~2 ( f 2  -  g2)_1[(l  -  f2)G3  +  (f4  -  g2)G]A3* 


+  (Ej_  +  E2G4)A2' 


Equations  (102)  and  (103)  together  with  (90),  (91)  and  (92)  yield 


(2)  «  %  [f_1(M2  +  2f2)A*e'f?  +  2gB*e-gc  +  E,e"Gc] 


(104) 


<rCn(2)  -  y  [2A*e‘fc  +  (1  +  g2)  B2*  e"9c  +  e"Gc] 


(105) 


—  (2)  _  w2  r_,-l 
nn 


y  C-f^d  +  g2)A*e-fc  -  2gB2e"9!:  +  Ege'G?] 


(106) 
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vTn'<4)  -  U2  *  t2)  7™'  ‘  *  ,  V  7^  .  r  7TT"  •  t  Ein  tit'-  0 


where  P  =  1  -  I2/J2 


(108) 


28 


The  complementary  solution  of  equation  (108)  is 

( 1 1 c  *  ★  *  * 

V  ■  a2coshjc  +  b2sinhjc  +  c2coshk?  +  d2  sinhkc  (109) 

Assume  the  particular  solution  of  equation  (108)  is 


V<1)P  =  A%~F<’  +  A!eF; 
4 


A5e 


(110) 


Substituting  (110)  into  (108)  we  can  obtain  A^  and  A*. 


I-  *9  *  E10FI  Ai 

4  '  F4  -  U2  I  k2)F2  7  A2 


★ 

A5  =3 


(F  %  -  F  E10)B1 
r  -  (j2  +  k2)F2  *  j2k2 


The  solution  of  equation  (108)  is 


( 1 )  *  *  *  * 

V  =  a2coshjc  +  b2sinhjc  +  c2coshkr,  +  d2sinhkc 


+  A*e~F!:  +  A%Fc 
4  5 


(111) 


Substituting  into  (107)  gives 


(]  )  .-1  .  *  * 
it  '  -  j  (a2sinhjc  +  bgcoshjg)  +  k(c2s1nhkc  +  d2  coshk^) 


+  E  e’Fs  +  E  eF^ 
L11  e  +  Li2  e 


(112) 


where  E 


F3  * 


..-2 


■lkF-.* 


*2 

1  b  Y1  r-2,' 


11  =  ~2T  a4  ■  (J'  +  p  )FA4  +  ~TTn  —f  rAl 


JP 


IT  P  Cr 
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I2/  C,2  1 


4  *;♦  u-*  ♦»“>«;♦ -a- 

J  p  I  J  P  C, 


*2 

i  6  ri  ■ 

77  Tf  1 


Equations  (111)  and  (112)  together  with  (85),  (86),  and  (87)  yield 

T“(1)  s  6r  C'  (M2j-1  +  2J2j) (a^sinhjc  +  b*coshjc) 

"  o 

-  2J2k(c*sinhkj;  +  d2coshk/;)  +  E^e  Fi*  +  E14eF^]  (113) 

a  ^  =  6J2  -5—  [2(a!coshjc  +  b*sinhjs)  +  (1  +  k2)(c?coshk^ 

5n  ‘  c 


*  ...  _  -Fr  _  Fr-i 


■15 


-  feu  -  a4 


■16 


FE12  +  A5 


■17 


*2 

9  *  9  9  b  Yi  * 

-  TFA4  -  11*  -  2J^)En - f-  Aj 

c2 


■18 


9  *  9  9  b  Yi  « 

i  FAg  -  (I  -  2^)E12 - f-  Bx 


Applying  the  boundary  conditions  and  the  continuity  conditions  at  the 

surface  layer/ substrate  interface,  the  six  constants  A2>  B2>  a2,  b*, 

★  * 

c2  ,  and  d2  can  be  solved  by  six  algebraic  equations 

2a*  +  (1  +  k2)c*  *  Eig 
j^d  +  k2)b*  +  2kd*  =  E20 

.-1  *  -1  *  *  *  -1  *  -fen  *  .ncn 

J  +  j  Cjb2  +  kskc2  +  kckd2  +  f  AA2e  s  +  gB2e  gsu  =  E21 


*  ■*  *  -*  *  -fsD  _  D*.-gsO 


Cja2  +  Sjb2  +  CkC2  +  Skd2  ’  A2e 


-  B2e 


"22 


2c j.a2  +  2s jb*  +  (1  +  k2)ckc2  +  (1  +  k2)skd*  -  — - ?  A* 

u  5J 

- 11  ±  gV930  B*  s  E 

t?  2  23 

j_1(i  +  k2)Sj.a2  +  j_1(l  +  k2)cJb*  +  2kskc*  +  2kckd* 


.  (1  +  g2)e-fsD  *  ,  2ge"gsD  *  _  _ 
^57 - A2+^?-B2-E24 

where  E^g  *  -  (E^  +  E^g) 
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*■ 


★ 


A 


* 


* 


* 


a 


where  a*  =  A. 

For  s  <  0  we  can  obtain  similar  solutions. 

*  *  *  *  *  ★ 

The  expressions  of  A^,  A2»  A3,  a4,  Ag,  and  Ag  are  given  in 
Appendix  II. 
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CHAPTER  III 


NUMERICAL  SOLUTION 


Previous  sections  have  shown  the  general  expressions  of  the 
stress  and  temperature  fields  in  the  transformed  space.  By  applying 
the  inverse  Fourier  transform  and  following  through  the  simplifica¬ 
tion  of  the  integrand,  we  can  obtain  results  for  several  cases  numeri¬ 
cally.  Here,  a  discussion  of  how  the  thickness  of  the  surface  layer 
affects  the  thermostress  and  what  material  parameters  govern  condi¬ 
tions  leading  to  crack  formation  will  also  be  addressed.  (Figures 
2-13). 


3.1  Mechanical  Stress  Field 


(1) 


EE 


*  Re  j(2ir ^ f  a^e"*S?  ds 


=  Re. 


C 5 ( M2 j  1  +  2J2j)  UjSinhjns  +  b*  coshes) 


+  2SJ2k(c^sinhkns  +  d^coshkns)]e”1s^  ds 

*  *  -isr  pe"(^+9)sD  (iy* 

(a1  sinhjns  +  b^oshjns)  e  ^  *  - - -  < — ^  (Bj^  +  B2CjCk 

+  B,c.s.  +  B.s.s.,  +  Bcc.s4)  — — *•  (B^c.s.  +  B,c  .c.,  +  B0c,s  . 

3jk  4jk  5kj  bjk  7jk  8kj 


+  BgSjSk)]  sinhjns  + 


iuf 

;^(cisj 


ck +  Wk 
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+  CgS^Cj)  coshjnsj  (coss^  -  isinss) 
|Re  (a*  sinhjns  +  t^coshjnsje  1S^| 


pe-(f+g)so  ( 

— r— | 


~~2  sinsc  (Bx  +  B2cjck  +  B3cjs|c  +  B4$jSk 


+  B5cks  .)  -  C0S-|£  (BgCjSk  +  B7CjCk  +  B qC k s j 

6  J 


+  Bg  Sj  sk ) I  si nh jns  + 


-  -77  s1nsE  (clsjck  *  c2sjsk 
5  J 


♦  *  c‘c<c'-1  C6  sjsk  *  C7sjck  *  C8CjCk 

6J 


"3cjsk  T  VjV 


*  Vkcj> 


coshjns 


)  5o"^f+9^sD  (Wf 

jTlSr  Sin$C  [(B1  +  B2C jck 


)  “  (sj 


+  B3cjsk  +  B4sjsk  +  B5ck sj) sinhjns  -  (c^  +  c 2sjS|c 
+  c3c  .sk  +  c4c.ck) coshjns]  -  -~|i  [(BgCjSk  +  B?cjc|c  +  BgCkSj 
+  BgSjsk)sinhjnS  ♦  (c5  +  CgS.sk  +  c7sjCk  ♦  c8c.ck 


+  cgskCj)coshjns]| 


a  £® - - - |_f_.  sins^tB. sinhjns  +  B»c.  (c. sinhjns  -  s  .coshjns) 

&  |5J2  1  ^  J 

+  B3sk(c jsinhjns  -  SjCOShjns)  +  S4sk(s jSinhjnS  -  cxo$hjns) 

+  Bcct [$. sinhjns  -  c. coshjns)]  -  cos|£  [Bfis.  (c  .sinhjns 
5  k  j  j  5J  J 

-  SjCoshjns)  +  B7ck(CjSinhjns  -  SjCoshjns)  +  BgCk(SjSinhjnS 

-  CjCoshjns)  +  Bgsk(s ^sinhjns  -  CjCoshjns)  -  c5coshjns]j 
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_  pe~^f+g^sD  |  Vf  sins?  [Bj^sinhjns  -  (B2ck  +  B3sk) si nh(D  -  n) Js 
A 

-  (B4sk  +  BgCk)cosh(D  -  n) js]  +  £££|i  L ( Bg sk  +  B7ck) 

5J 

si nh (D  -  n) Js  +  (BgCk  +  Bgsk)  cosh(D  -  n)js  -  c5coshjns]j 

*  *  -is?  pe-(f+g^sDir  i 

(CjSinhkns  +  d^oshkns)  e  *  * - - -  - 2  ^diskcj 


*  d2‘j‘k  *  d3*j‘k  +  VjV  *  ^  (V,*l  *  Vj‘k 


*  VA  *  d8skc  j  *  V  si"hk',s  + 


iuf 


ck  +  P2SjSk  +  P3skcj 


*  P4C4Ck>  *  ^  ‘Vd'k  *  '.Vi  *  VA  *  Vj‘k  *  V 

coshkns j  (cossc  -  isirts?) 


Re  ac^slnhkns  +  dj  coshkns)  e 


Pe~(f+g)sD  j  coss5 


S  j  77*  [(VkC3  +  d2CJCk  +  d3S3Ck  *  VoV 

sinhkns  -  (P5c..ck  +  PgS^.  +  PyS^  +  PgS^  +  Pg)coshkns] 


+  — 2  sins?  C(d5SjSk  +  d6SjCk  +  d7cjc|C  +  d8skcj  +  <V3inhkr’s 

6J 

+  ^Plsjck  +  P2sjsk  +  P3skcj  +  P4C jck ^coshkr*S^  | 

5  -(f+g)sO  {  _ 

»  ~ - \  -co--^  +  d4Sj)cosh(D-ri)ks 


uf  .... 


+  (d0c.  +  d,s.)sinh(D  -  n)ks  +  P0coshkns] - *•  sins£[(d,.s. 

t  J  J  J  3  x.l6  3  J 
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+  dQc.)  cosh(D  -  n)ks  +  (d,-s.  +  d,c.)sinh(D  -  n)ks 
a  j  o  j  /  j 

-  d„sinhkns] 


\ 


5  e"(f+9)sD  2-1  2 

- -  S  [ 6(M  j  +  2JZj) 


Therefore  -J-  -K  f 

55  /2ir  Ko  sJ^  yQ 

jufsins£  [  -  B^sinhjns  +  (B2ck  +  B3sk)sinh(D  -  n)js  +  (B4sk 

+  BgCk)cosh{0  -  n)js]  -  coss£[(BgSk  +  B7ck)sinh(D  -  n)js 

)  2  ( 

+  (BgCk  +  BgSk)cosh{D  -  n)js  -  CgCoshjns]j  +  26 Jk  jcosss 
[(d^Cj  +  d4Sj)cosh(0  -  n)ks  +  (d2Cj  +  d3Sj)sinh(D  -  n)ks 
+  Pgcoshkns]  +  u^sinsc  [(dgS^  +  dgCj)cosh(D  -  n)ks 


+  (dgSj  +  d7Cj)sinh(D  -  n)ks  -  dgsinhkns]  j 


ds. 


Similarly, 


(1) 

\r\ 


.  2  u2  / 
Po  1 


P  e 


■(f+g)sD 


2| y^coss?[B jCoshjns  +  ^B2ck  +  B3sk^ 


cosh(D  -  r») js  +  (B4sk  +  BgCk)sinh(D  -  ri)js]  +  sinss[(Bg$k  +  B7ck) 
cosh(0  -  n)js  +  (Bgck  +  Bgsk)sinh(D  -  n)js  +  c5sinhjns]j 

+  (1  +  k2)  j-  sinscC(d^Cj  +  d4Sj)sinh(0  -  n)ks 
+  (d2Cj  +  dgSj)cosh(D  -  n)ks  -  Pgsinhkns]  +  ufcoss£ 

[(dcs.  +  dQc.)sinh(D  -  n)ks  +  (d,s.  +  d,c.)cosh(D  -  n)ks 

DJOJ  oj/j 


+  dgcoshkns] 


ds 
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(1)  _  2  y2  /  P  e~(f+9)sD 


=  J_^2  f 

m  PoJ 


ijj_1{l+k2) jufsinss  [B^inhjns 


-  (B2ck  +  B3sk)sinh(D  -  n)js  -  (B4sk  +  B5ck)cosh(D  -  n)js] 

+  cossc[(BgSk  +  B7ck)sinh(D  -  n)js  +  (BgCk  +  B^) 
cosh(D  -  n)js  -  c5coshjns]|  +  2k  j  -  cossc  C(djCj  +  d 4s ^ ) 
cosh(D  -  n)ks  +  (dgCj  +  dgSj)  sinh(D  -  n)ks  -  Pgcoshkns] 

-  sinss  HdgSj  +  dgCj)cosh(D  -  n)ks  +  (d^s^  +  d7Cj)  sinh(D  -  n)ks 

+  dgSlnhkns]  j  j  ds 

U)  ,  J_  “2  1  f  £  L  f-l(H2  „  2t2)e-'3^  *  fc> 

^  ^  5?  J0  A  I 

[yfsinsc(H1Sj  +  HgCj  +  HgCk  +  H4sk)  +  coss^HgC^  +  HgS^ 

+  H?sk  +  Hack)]  -  2g  e'(fsD+9c)CUfsinsc(L1sj  +  LgC^  +  L3ck 
t  L4skl  +  cossetLsCj  ♦  L6Sj  *  LySk  *  L3ck)]  j  ds 

.<«.  ±  £  -ij  f  I  S  |  Z'-W  *  f5>  [„fcosst  Wl. .  *  HfC, 

{n  /sr  po  «o2  y0  1  i  f  i  d  2  j 

*  H3ck  *  H4Sk’  •  sios*lH5Cj  *  H6sj  *  H7sk  *  Vk'1 

-  (1  +  g2)  e"(fsD+9?)  Cufcoss?(L1Sj  +  LgC^  +  L3ck  +  L4sR ) 

-  sinss  (L5Cj  *  LgSj  *  L,sk  *  L,ck)]J  ds 
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£  1  fi  s  J  f-Ml  *  92)  e-(9sD+f‘> 
nn  m  po  jj2  J0  a  t 

[ufsinsc(H1Sj  +  H2Cj  +  H3ck  +  H4$k)  +  coss^HgCj  +  HgSj  +  H7sk 

+  Hgck)3  -  2ge"(fsD  +  9?)  [ufstnsc(L1sj  +  Lzc.  +  L3ck  +  L4sk) 

+  coss?  (LgC.  +  LgSj  +  L?sk  +  LgCk)]  |  ds 

The  above  integrals  are  too  complicated  to  integrate  directly; 
therefore  a  numerical  integration  technique  will  be  employed.  There 
are  many  numerical  schemes  for  integration.  Here  subroutine  QUANC8 
will  be  used  which  is  based  on  the  8-panel  Newton-Cotes  rule  [8].  In 
the  integrals,  s  «  0  is  a  singular  point.  The  method  used  to  avoid 
the  singularity  is  discussed  in  Appendix  III. 

3.2  Temperature  Field 


(F  +  0G)e"FnS+  (F  -  BG)e(rr2D)Fs 
F2(l  -  e~2,fsD)  +  gFGU  +  e”2FsD) 
1/2 

Let  F  *  (1  -  iRj^/s)  =  Y^cosej^  +  isin^) 

1/2 

G  *  (1  -  iR2/s)  =  Y2(cose2  +  isine2) 
where 

yl  »  (1  +  Rx2/s2) 1/4,  9j  -  1/2  tan_1(-R1/s) 
y2  3  (1  +  R22/s2)1/4,  02  -  1/2  tan_1(-R2/s) 
F2  (1  -  e"2FsD)  +  bFG(1  +  e‘2FsD)  -  Tj  +  iT2 


r  — ★ 

7^  =  2- 

v  s 


where 


T 


1 


cos  29^ 


-2y, sDcoso, 

1  -  cos  (2Y1sDsin91)  e  1 
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-2Y.SDCOS9,  )  / 

-  si n291  si n(2Y2 sDsi nex )  e  +  By^  )  cos( e1  +  02) 


( 


1  +  cos  (2Y1sDsine1)  e 

“2YJSDCOS0,) 

6  1 


-2y1sDcos0^ 


+  si n ( 0^  +  92 ) s i n  (2Y1sDs1n91) 


2  l 

r2  -  Yi  !  Sin 20, 


1  -  cos  (2Y1sDs1n01)  e 


-2y1sDcos9^ 


+  COS20 


1 


-2y1sDcos01|  j 

sin(2Y1sDsine1)e  +  8y1y2  sin(91  +  ^  ^  +  cos(2Y1sDsine1) 


-2YTSDCOS0,  ) 

g^Yj^sDcosBj^]  -  cos(91  +  a2)  sin(2Y1sDsin91)  e  1 
(F  +  8G)e“FnS  =  T3  +  iT4 


where 


-YiSncose. 

T3  =  e  [cos(Y1snsin01)(Y1coso1  +  8y2cos92)  +  s^W^nsin©^ 

(Y1sine1  +  0Y2sin02)] 

-Y^sncose^ 

T4  =  e  1  [cos(Y1snsinei){Y1sine1  +  BY2sine2)  "  sin^snsingj) 

(y1cos01  +  By2cos02)] 

(F  -  bG  )  e(n"2D)Fs  =  T5  +  lTg 
where 

i  r2D)sY1cos8,  j 

T5  «  e  cos[(n  -  2D)SY1sin91](Y1cos81  -  by2cos92) 

-  sin[(n  -  2D)SY1sin913(Y1sin91  - 
(WObY^ose.  | 

T{  s  e  j  cosCU  ~  2D)SY1sin91](Ylsin91  -  8Y2sine2) 


8Y2sin02) | 
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+  si n[ ( n  -  2D) s^sinejKYjCos^  -  by2C0S92^ 


I 


,(1)  _  0. 

*  s 


(T3  *  tT4)  *  <TS  ♦  iT6) 


T-nr2 


f  (T7  *  V 


where 


T^Ij  +  Tg )  *  T2(T4  *  T6) 

T7  t"5  77  2 

'l  ‘2 


T8 


T1(T4  +  T6)  “  T2(T3  +  V 

T  2  +  t  2 
‘l  2 


M)  = 


Re  (J_  f  ;<»»  e'is5  ds) 
)  /2ir  -/  I 


-  Re  | 

00 

^Jo 

=  _2_ 

CO 

(’■ 

Similarly, 

.(2)  .  2 

♦  ~ 

00 

I5*' 

where 

T  _  T1T9 

+  T2T10 

+  iTg}  e~is^  ds 


pcosss  +  TgSins^)  ds. 


ir 


12* 


ii  ~JtT7t 

*1  *2 


.  _  T1T10  "  T2T9 

12  i  2  +  T  2 
'l  1 2 
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T9  =  2V 


C-Y1sDcos91+(D-n)sY2cos02] 


|  cos91  cos  [-Y^sDsine^ 


+  {D  -  n)sY2sin02]  “  si n 91  si n  Or^sDsin©!  +  (D  -  n)sY2sin023  j 
C-y1sDcos01+(D-t1}sy2cos92]  ( 

Tw  *  2yj  e  jcosOjSin  [-YjSDsinOj 

+  (D  -  n)sY2sin©2]  +  sinOjCos  [-Y^DsinsOj^  +  (D  -  n)sY2sine23 
3.3  Thermal  Stress  Field 


,(1)  - 

55 


■i?  ■/>■* 


where 


F]L  =  -  (M2j-1  +  2J*j)  [cossc  |  -EjJ  C-AjSinhjns  +  (A3ck  +  A5sk) 
si nh(D  -  n)js  +  (A7sk  +  Agck)cosh(D  -  n)js]  -  (1  +  k2) 

[(BgSk  +  B^ck )  sinh(D  -  n)js  +  (B^cfc  +  B^sk)cosh(D  -  n)js 

-  (B10Sk  +  BllCk  +  B12sj  +  B13Cj)  S1'nhjnS  +  (B12C j  +  B13s j} 
coshjns]  +  2(cJ  +  cjck  +  cjsk Jcoshjns j  +  sins?  j  -  E|g 
[-  AjSinhjns  +  (A3ck  +  Ag sk ) si nh (D  -  n)js  +  ( A? sk  +  Agck) 

cosh(D  -  n)js]  -  (1  +  k2)[(BgSk  +  B^c k) si nh(D  -  n)js  +  (BgCk  +  Bgsk) 

cosh(D  -  n) js  -  (B11Qsk  +  B^ck  +  b|2Sj  +  B^3c j )  sinhjns 

+  (812cj  *  6l3>  coshjns]  +  2(c|  +  c2ck  +  c3sk)  coshjnsj] 


i2  •< 


I- 
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+  2J2k  [cossC|2[(B^0S|c  +  Bj1ck)sinhksn  +  +  djs^}  cosh(D  -  n)ks 

+  (djcj  +  dgSj jsinh  (D  -  n)ks  +  (B^Sj  +  B^Cj)  sinhkns  +  (BjQck 

+  B]1sk)coshkns]  +  Ejg  [(AgSj  +  AgC^cosh  (D  -  n)ks  +  (A10s..  +  A^.) 

si nh ( D  -  n)ks  -  A2sinhksn]  +  (1  +  k2)(Pjcj  +  Pjsj  +  Pj)  coshkns| 

+  sinsej2[(B110sk  +  B’1ck)  sinhkns  +  (d^c ^  +  dgS^)  cosh(D  -  n)ks 

+  (dj  +  dgSj)  si  nh  ( D  -  n)ks  +  (B^Sj  +  BjgCjJsinhkns  +  (B|Qck 

+  Bjxsk)  coshkns]  +  Ejg  [(AgSj  +  AgCj)cosh(D  -  n)ks  +  (A1QSj  +  A^) 

sinh  (D  -  n)ks  -  A2s1nhkns]  +  (1  +  )(p7cj  +  P8  sj  +  Pg)coshknsJ  ] 

+  T13cossc  “  T14sins? 

where  the  superscript  y  designates  the  real  part  of  that  constant,  i 
the  imaginary  part.  For  example,  E2Q  ■  E^  +  iE2Q  . 


T13  *  e 


*  [E^  cos  (r^nsine^  +  E13  sin  (yjSnsinej)] 


T  3  6 
14  6 


+  eTl  1  [Ej  cos  (y^snsinej)  -  E14  sin  (y-jSnsi ne^^ ) ] 
-v. sncose,  j 

i[E13  cos  (yjSnsinflj)  -  Ej  sin  (YjSnsinQj)] 


+  eYlSn  01  CeJ  cos  (y^sin^)  +  Ej  sin  (y^nsinOj)] 


,(D  -  -f^ds 

5n  Jn  *  2 


43 


wfiere 


F2  =  2Csinsc  |  E19C(A3ck  +  A5sk)cosh(D  -  n)js  +  (A?sk  +  AgCk) 

o 

s1nh(D  -  n)js  +  A^oshjns]  +  (1  +  k  )[Bgysk  +  B7yck) 
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CHAPTER  IV 


CONCLUSION 

In  the  rapid  traversing  friction  source  model,  the  thermal  effect 
dominates  the  stress  field  that  eventually  leads  to  failure.  Numerical 
results  are  presented  in  the  text  corresponding  to  five  different  cases 
of  material  properties  and  geometry  to  study  the  layer  thickness  effect, 
the  stiffness  ratio  effect,  the  thermal  diffusivity  ratio  effect,  and  the 
effect  of  pressure  distribution.  Only  the  stresses  at  £  >  0  are 
discussed,  since  the  maximum  stress  occurs  at  0  <  £  <  1  in  all  the 
cases.  Case  1  (Figures  2-4)  illustrates  the  effect  of  a  soft  layer  for 
which  the  material  of  the  surface  layer  is  represented  by  zirconium;  the 
material  of  the  substrate  is  stellite  III  which  is  harder  than  zirconium; 
the  dimensionless  layer  thickness  D  is  two;  and  the  pressure  distribution 
is  uniform.  Figure  4  shows  that  the  maximum  thermal  stress  for  case  1  is 
1.1  P0  .  Case  2  (Figures  5-6)  illustrates  the  thickness  effect,  for  which 
the  materials  of  the  surface  layer  and  the  substrate  are  the  same  as  in 
case  1,  but  D  is  different.  Figure  6  shows  that  the  maximum  thermal 
stresses  are  1.25  p0  ,  1.6  p0  ,  and  1.05  P0  for  D=  0.01,  0.1,  and  2.0, 
respectively.  The  phenomenon  indicates  that  the  stress  is  maximum  when 
the  layer  thickness  is  in  the  neighborhood  of  *7=  0.1,  the  depth  of 
maximum  thermal  gradient.  Case  3  (Figures  7-8)  illustrates  the  effect  of 
thermal  diffusivity  ratio,  in  which  the  material  properties  are  the  same 
as  those  in  case  1  except  that  the  thermal  diffusivity  of  the  layer 
material  is  one-half  that  of  zirconium,  with  D=2.  In  this  case,  both  the 
temperature  and  the  thermal  stress  decrease  in  comparison  to  those  of 
case  1.  In  case  4  (Figures  9-10)  the  effect  of  a  hard  coating  is 
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illustrated,  in  which  the  material  properties  are  the  same  as  those  in 
case  1  except  that  the  Young's  modulus  of  the  layer  material  is  five  times 
that  of  zirconium,  with  D=  2.  It  is  shown  in  Figure  10  that  the  thermal 
stress  is  much  larger  than  that  in  case  1.  In  case  5  the  effect  of 
pressure  profile  is  studied.  In  this  case  (Figures  11-13)  the  pressure 
distribution  is  parabolic  and  the  total  pressure  is  the  same  as  the 
uniform  pressure  case  and  the  material  properties  are  the  same  as  in 
case  1.  The  thermal  stress  due  to  this  type  of  excitation  is  larger  than 
that  with  uniformly  distributed  pressure. 

Base  on  the  cases  studied,  it  may  thus  be  concluded  that  a  stiff  surface 
layer,  which  is  less  compliant,  would  result  in  higher  thermal  stress. 
Surface  layer  materials  with  low  thermal  diffusivity  are  generally  a  result 
of  high  thermal  capacity  which  in  turn,  results  in  less  temperature  rise, 
thus  lower  thermal  stress.  Finally,  the  effect  of  surface  layer  thickness 
depends  on  the  thermal  layer  which  is  defined  by  the  depth  of  maximum 
temperature  gradient.  Therefore,  substantial  reduction  in  thermal  stress 
can  be  achieved  by  decreasing  the  thermal  diffusivity,  by  decreasing  the 
modulus  of  elasticity  of  the  surface  layer,  or  by  designing  the  layer 
thickness  such  that  the  interface  is  away  from  the  depth  of  maximum 
temperature  gradient,  which  in  the  present  numerical  problem  is  77=0.1. 
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Figure  2.  Dimensionless  Mechanical  Principal  Stress 

The  material  of  the  surface  layer  is  zirconium.  The 
material  of  the  substrate  is  stellite  III.  D  =  tt/l 


Thermal  principal  stress/average  pressure 

Combined  thermomechanical  orlnclpal 
stress/average  pressure  (uniform  pressure) 
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The  materials  of  the  surface  layer  and  the  substrate  are  the  same 
as  Figure  1.  D=2,  1=0.01  in. 
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substrat 


Combined  thermomechanical  principal  stress/ 
average  pressure  (uni  lorn  pressure) 
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Dimensionless  Mechanical  Principal  Stress 


too 


The  materials  of  the  surface  layer  and  the  substrate  are  the  same  as 
Figure  1.  D  =  H/£  =  2,  £  =  0.01  in. 


Combined  thermomechanicat  principal 
stress/average  pressure  (parabolic 


The  materials  of  the  surface  layer  and  the  substrate  are  the  same 
as  Figure  1.  D  =  H/l  =  2,  t  =  0.01  in. 
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P1  *  ”A8*  P2  *  'A10*  P3  =  ”A4’  P4  =  ”A6*  P5  =  A3/k*  P6  =  A5/k* 


Py  =  A7/k,  Pg  =  Ag/k,  Pg  =  Ax/k 
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2kg  -  2(1  +  k  ) 


+  k2)  +2U_Mcf).+ 


a  aii-t  k2i  -  2(-  +  (1  +  k2)2  +  -  2(1  +  k2) 

J  6J  dJ 

.  k(l  +  k2)  (1  +  k2)(l  +_  g2)  _  2(1  ^  g2)  .  2(1  +  k2}  +  4 

J  Zi2  x/i2 


(1+k2) 


0  j”1 ( 1+k2)  0 


j_1s.  rXc.  ks. 

J  J  N 


sj  ck 


2Sj  (l+k2)ck  (l+k2)sk 


j"1(l+k2)$j  j"1(l+k2)cj  2ksk  2kck 


(l+g2)e-fsD 


e"fsD  — (L,  s  .  +  L,c.  +  L,c.  +  L.sJ  -  (Uc,  +  L,s. 
6J2  1  J  i  J  J  K  5>T  3  J  J 


+  L7  Sk  +  L8Ck) 


Lx  »  -  2(1  +  k2)  +  -^  +  (1  +  k2)  -  2(1-*g-  ■ 


.  ,  -k(l  +  k2)  f_  2(1  +  k2)  .  (1  +  k2)(l  +  g2)  .  4  _  2(1  +  g2) 
2  J  f  fdJ2  7  f  6J2 
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APPENDIX  II 


THE  EXPRESSIONS  OF  A*,  Ag,  A*,  aJ,  A*,  AND  Ag* 


E19 

0 

(1+k2) 

0 

0 

0 

E20 

j'Ml+k2) 

0 

2k 

0 

0 

E2I 

fS 

ksk 

kck 

f-Vfs0 

-gsD 

ge 

* 

Al  = 

E22 

sj 

ck 

sk 

-fsD 

-e 

-e‘9S0 

2sj 

(l+k2)ck 

(l+k2)sk 

-2e-fsD 

-(l+q2)e"9sD 

E23 

“T7~ 

SJ2 

E24 

j"l(l+K*)Cj 

to 

CM 

2kcic 

(U»2)e-fsD 

f5J2 

2(je"9sD 

6J2 

*  e-(f*9ls0  tEu(8j  *  BjCjC*  *  Bj^Cj  *  *  BjSjC,.) 

♦  »  ♦  k2HB6skCj  ♦  B7CjCk  ♦  BgS jCk  *  BgSjsk  ♦  B10sk  t  Buck 

+  B12Sj  +  ®i3c j)J 


where 
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E23  >1+k  lck  ll+k  )sk 


2cj 


- -y~  -  a  J 

1  C.  K 


4J 


j*l(l+k2)s. 


E24  2ksk 


2kc„ 


’k  ffij2" 


6J 


r 


-lf+g)sD 
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+  c7cjck)  +  {1  +  k  ,(c8sjsk  +  C9S jck  +  c10cjck  +  cllskcj 


+  C12Cj  +  c13s j  J 
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E19 
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E20 
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+  dllcjck 
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where 
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4 


d 
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2k 
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13 

10 

4 

6 

2 


78 


2 

0 

(1+k2) 

E19 

0 

0 

0 

j*1(l+k2) 

0 

E20 

0 

0 

II 

*<r 

,  J-S 
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APPENDIX  III 


THE  METHOD  USED  TO  AVOID  THE  SINGULARITY  IN  THE  INTEGRAL 
In  this  appendix,  the  method  used  to  avoid  the  numerical  singu¬ 
larity  of  the  integrals  encountered  in  the  analysis  is  discussed  by 

f4  1 

means  of  the  example  integral  /  —  dx.  This  integral  has  a  numerical 

J0  /x 

singularity  at  x  =  0,  and  therefore  cannot  be  integrated  numerically 
from  zero  to  four.  One  approach  to  avoid  the  singularity  is  to 
replace  the  lower  limit  of  the  integral  by  e  which  is  an  arbitrary 
small  number  greater  than  zero.  Integration  is  carried  out  with 
decreasing  values  of  e  until  a  value  of  the  integral  is  obtained  which 
is  considerably  different  from  the  previous  integral  value  correspond¬ 
ing  to  the  previous  value  of  e.  The  last  value  of  the  integral  thus 
obtained  is  taken  to  be  the  best  approximation  to  the  actual  value  of 
the  integral.  For  the  example 


1_ 

/x 


dx 


the  exact  solution  is  4.  The  results  of  the  numerical  computation  of 
A  for  various  values  of  e  are  tabulated  below. 
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e 

A 

10‘9 

3.99994 

o 

1 

o 

r— • 

3.99998 

io~n 

4.00000 

io"12 

4.00005 

10-13 

4.00018 

io'14 

4.00063 

io“15 

4.00204 

io-16 

4.00647 

T-l 

1 

o 

4.02051 

10~18 

4.06490 

IO'19 

4.20526 

10 -20 

4.64914 

The  table  indicates  that  the  numerical  integration  diverges  for  values 
of  e  smaller  than  IO"*8.  e  *  10"18  therefore  gives  the  best 
numerical  solution  to  the  problem. 
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